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Abstract. Anomalous X-ray Pulsars (AXPs) belong to a class of neutron stars believed to harbor the strongest magnetic fields 
in the universe, as indicated by their energetic bursts and their rapid spindowns. We have developed a theoretical model that 
takes into account processes in the atmospheres and magnetospheres of ultramagnetic neutron stars, as well as the effects of 
their strong gravitational fields on the observable properties. Using this model, we have analyzed the X-ray spectra of a number 
of AXPs. We find that in all cases, the X-ray spectra are described very well with this emission model. The spectroscopically 
measured magnetic field strengths of these sources are in close agreement with the values inferred from their spindown 
properties and provide independent evidence for their magnetar nature. The analysis of spectral data using this physical model 
also sheds light on the long-term evolution of AXPs. 
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INTRODUCTION 

Anomalous X-ray Pulsars (AXPs) are thought to be the 
observational manifestations of a class of ultramagnetic 
(B> 10 14 G) neutron stars, also called magnetars (see 
Woods & Thompson 2006 and Kaspi 2006 for reviews 
on magnetars and AXPs). Among the numerous spectral 
and timing properties of these isolated X-ray sources, 
two stand out for our focus in these proceedings. Their 
X-ray spectra are soft but non-Planckian, traditionally 
described by empirical functions such as a blackbody 
(kT ~ 0.3—0.6 keV) plus a power law (with photon 
index F ^2.5—4) and, less frequently, by a sum of two 
blackbody functions (see, e.g., Gotthelf & Halpern 2005; 
Kaspi 2006). The second property is their high spin- 
down rates, with P ~ 10~ n s s . 

A convincing, albeit indirect, argument for their strong 
magnetic fields arises from these large spindown ratesj 
(e.g., Kouveliotou et al. 1998). Assuming the neutron 
stars spin down due to magnetic braking of a dipole in 
vacuum, their magnetic field strengths can be estimated 
by B dip = 2.48 x 10 14 (P/6 s) l / 2 (P/10- n s s" 1 ) 1 / 2 G, 
for a neutron star moment-of-inertia / = 10 45 g cm 2 and 
a neutron star radius of R = 10 km. The dipole fields 
associated with AXPs thus exceed 8>5x 10 13 G. The 
dipole spindown formula makes numerous assumptions 
when connecting period derivatives with a magnetic field 



The energetics and the timescales of intense, super-Eddington, ran- 
dom bursts of X-rays or soft gamma-rays seen in AXPs and the closely 
related Soft Gamma-ray Repeaters that last a fraction of a second 
also suggest independently the existence of very strong magnetic fields 
(Thompson & Duncan 1995). 



strength, such as a fiducial angle between the magnetic 
and rotation axes and the absence of other torques on 
the neutron star (Spitkovsky 2006). The dipole magnetic 
field inferred in this way has never been compared with 
an independent, spectroscopic measurement for an iso- 
lated pulsar. 

Recently, there has been significant theoretical work 
on the emission from the atmospheres and magneto- 
spheres of magnetars. As part of these efforts, we have 
developed a physical model of emission from a magne- 
tar that takes into account processes in its atmosphere as 
well as in its magnetosphere. The Surface Thermal Emis- 
sion and Magnetospheric Scattering (STEMS) model is 
based on the radiative equilibrium atmosphere calcula- 
tions presented in Ozel (2003) but also includes the ef- 
fects of magnetospheric scattering of the surface radi- 
ation as discussed in Lyutikov & Gavriil (2006) and 
Giiver, Ozel & Lyutikov (2006). We also take into ac- 
count the general relativistic effects in the strong gravi- 
tational field of the neutron star, making our models di- 
rectly comparable to the wealth of spectral and timing 
data on AXPs. Naturally, comparison with such data is 
the ultimate test of any theoretical model. At the same 
time, a model that can describe consistently and in de- 
tail the spectra of AXPs can be used to understand the 
physical properties of these sources and their emission 
mechanisms. 

In this proceedings paper, we present the results of ap- 
plying the Surface Magnetospheric Scattering and Sur- 
face Emission Model to the soft X-ray data of four AXPs. 
In particular, we measure the magnetic field strength 
of these sources spectroscopically and we investigate 
the connection between the spectroscopically determined 



magnetic field strengths with those inferred from dipole 
spindown. 

THE THEORETICAL MODEL 

In highly magnetic, ionized neutron star atmospheres, 
polarization-mode dependent transport of radiation that 
includes absorption, emission, and scattering processes 
determines the continuum spectrum (see, e.g., Ozel 
2001, 2003). Furthermore, the interaction of the pho- 
tons with the protons in the plasma gives rise to an ab- 
sorption feature at the proton cyclotron energy E p = 
6.3 (B/10 15 G) keV. This absorption feature is weakened 
by the vacuum polarization resonance, which also leads 
to an enhanced conversion between photons of different 
polarization modes as they propagate through the atmo- 
sphere. 

In the magnetospheres of magnetars, currents sup- 
porting the ultrastrong magnetic fields can lead to en- 
hanced charge densities (Thompson, Lyutikov, & Kulka- 
rni 2002), which reprocess the surface radiation through 
resonant cyclotron scattering (Lyutikov & Gavriil 2006; 
Giiver, Ozel, & Lyutikov 2007). We calculate this effect 
using the Green's function approach described in Lyu- 
tikov & Gavriil (2006) assuming that the magnetosphere 
is spherically symmetric and the field strength follows a 
1/r 3 dependence. 

In our spectral models, we include the relevant pro- 
cesses that take place on the magnetar surface and its 
magnetosphere, which depend only on four physical pa- 
rameters. The first two parameters, the surface magnetic 
field strength B and temperature T, describe the condi- 
tions found on the neutron star surface. The third param- 
eter denotes the average energy of the charges j3 = v e /c 
in the magnetosphere, while the last parameter is related 
to the density N e of such charges and indicates the opti- 
cal depth to resonant scattering by z = ajN e dz. Here, 
a is the cross-section for resonant cyclotron scattering. 
We also assume a fixed value for the gravitational accel- 
eration on the neutron star surface of 1.9 x 10 14 cm s~ 2 , 
obtained for reasonable values of the neutron star mass 
and radius. 

We calculated model X-ray spectra (in the 0.05 - 9.8 
keV range) by varying model parameters in suitable 
ranges that are in line with the physical processes we in- 
corporated into the models: surface temperature T = 0.1 
to 0.6 keV, magnetic field B = 5x 10 13 to 3 x 10 15 G, 
electron velocity j3 =0.1 to 0.5, and optical depth in 
the magnetosphere T = 1 to 10. From the set of calcu- 
lated spectra, we created a table model which we use 
within the X-ray spectral analysis package XSPEC (Ar- 
naud 1996) to model the X-ray spectra of AXPs. 

Our models predict strong deviations from a Planckian 



spectrum, with a hard excess that depends on the surface 
temperature as well as the magnetic field strength, and 
weak absorption lines due to the proton cyclotron reso- 
nance. Both the atmospheric processes and the magne- 
tospheric scattering play a role in forming these spectral 
features and especially in reducing the equivalent widths 
of the cyclotron lines. 

ANALYSES OF AXP SPECTRA 

In this proceedings paper, we present the analysis of a to- 
tal of four XMM-Newton observations of four AXPs. For 
4U 0142+61, 1RXS J1708-4009, and XTE J1810 197 
we chose the longest available X-ray observation carried 
out by XMM or Chandra observatories (i.e., the observa- 
tion with the highest total counts). For IE 1048.1-5937, 
we used the longest observation of this source in quies- 
cence. A longer observation taken just after a burst from 
this source will be presented elsewhere. 

In Table [T] we present the list of the archival pointed 
X-ray observations of each source analyzed in this study. 
All of these observations were taken with the European 
Photon Imaging Camera (EPIC) PN camera. The obser- 
vations of IE 1048.1-5937, 1RXS J1708-4009, and 
XTE Jl 8 10 - 197 were taken in the Small Window Mode, 
while the observation of 4U 0142+61 was taken in the 
Fast Timing Mode. 

The spectral analysis was performed using the XSPEC 
11.3.2.t (Arnaud 1996). We assumed a fiducial gravita- 
tional redshift correction of 0.2, which corresponds to a 
neutron star with mass 1 .4 M and R = 13.8 km. We cal- 
culate the fluxes for the 0.5 — 8.0 keV energy range and 
quote errors for 90% confidence level. 



AXP 4U 0142+61 

4U 0142+61 is the brightest known AXP and has 
historically been very stable. Rotating with a 8.69 s 
period (Israel et al. 1994), it spins down at a rate of 
P « 0. 196 s s" 1 , yielding a B dip = 1 .3 x 10 14 G using the 
dipole spindown formula. Multiple X-ray observations of 
the source showed a long epoch of nearly constant flux 
levels as well as a relatively hard X-ray spectrum (Juett 
et al. 2002; Patel et al. 2003; Gohler, Wilms & Staubert 
2005). Recently, the source exhibited SGR like bursts 
(Kaspi, Dib & Gavriil 2006; Dib et al. 2006; Gavriil et 
al. 2007) for the first time. 

4U 0142+61 has also been detected in hard X-rays 
with INTEGRAL (Kuiper et al. 2006, den Hartog et 
al. 2007a). The hard X-ray spectral component in the 
20 — 230 keV energy range is well described by a power 
law model of index 0.79 and the corresponding flux 



TABLE 1. Observations used for this study. 



Exposure 

Source Satellite Detector Mode Time (ks) Obs ID Obs Date 



4U 0142+61 XMM-Newton EPIC-PN Fast Timing 21.1 0206670101 Jul 25 2004 

IE 1048.1-5937 XMM-Newton EPIC-PN Small Window 32.44 0307410201 Jun 16 2005 

1RXS J1708-4009 XMM-Newton EPIC-PN Small Window 44.9 0148690101 Aug 29 2003 

XTEJ1810-197 XMM-Newton EPIC-PN Small Window 42.2 0301270501 Mar 18 2005 



is 1.7 x 1(T 10 erg cm^s" 1 (den Hartog et al. 2007a), 
which exceeds by a factor of ~2 the unabsorbed 2-10 
keV flux. If this component extends without a break 
towards lower photon energies it contributes significantly 
to the soft X-ray flux in the 7-10 keV range. Because 
of this, in our present analysis, we take into account the 
effect of this component by using the fits to the hard X- 
ray observations reported by den Hartog et al. (2007a), 
assuming that this component extends to the soft X-rays 
without a break. 



IE 1048.1-5937 

Several properties distinguish IE 1048.1—5937 from 
the other AXPs. An ongoing RXTE monitoring cam- 
paign (Gavriil and Kaspi 2004) revealed that it shows 
long-lived pulsed flux flares in addition to SGR-like 
bursts. The spindown of the 6.452 s period pulsar is very 
unstable, with period derivative values in the range P = 
0.8546(50) -3.81 x lO^'ss" 1 (Kaspi et al. 2001). This 
yields a large range of dipole magnetic field strengths es- 
timated from spindown. As a conservative range, we will 
adopt B dip = 2.4 - 4 x 10 14 G for this source. 

IE 1048.1—5937 has been observed by the Chandra 
and XMM observatories as part of ongoing campaigns to 
monitor the variability of this AXP. The longest observa- 
tion to date was taken by XMM on 16 June 2003, shortly 
after bursting activity. To focus on the quiescent proper- 
ties of this source, as with the other AXPs, we analyze 
here the longest observation in quiescence, taken on 16 
June 2005. 



1RXS J1708-4009 

1RXS J 1708 -4009 is an 1 1.0 s AXP, initially thought 
to be a fairly stable rotator (Israel et al. 1999). In 
the last several years, the source experienced multiple 
glitches (e.g., Dib et al. 2007) that interrupted stretches 
of steady spin-down. A period derivative of P « 1.4 — 
1.9 x 10~ n s s _1 yields a dipole magnetic field strength 
of B dip =4.0-4.7 x 10 14 G. 

As in the case of 4U 0142+61, 1RXS J 1708 -4009 ex- 
hibits a hard, pulsed hard X-ray tail extending to energies 



up to ~150 keV (Kuiper et al. 2006). Here we adopt the 
values given by den Hartog et al. (2007b) with F = 1.17 
and the 20-250 keV flux 6.2 x 10" 11 erg cm^s" 1 . 



XTE J1810-197 

In the opposite extreme from 4U 0142+61, 
XTE J1810 197 is the most variable confirmed AXP 
observed to date. It was discovered (Ibrahim et al. 

2004) in 2003 when it suddenly brightened to more 
than 100 times its quiescent value (Halpern & Gotthelf 

2005) during an outburst. The source showed a steady 
decline of its X-ray flux thereafter, down to unusually 
low quiescent flux levels that have been determined 
from archival XTE and ROSAT data, accompanied by 
significant spectral changes (Gotthelf & Halpern 2006), 
earning it the title of the transient AXP. The detection of 
characteristic X-ray bursts (Woods et al. 2005), similar 
to those seen in other AXPs (Gavriil, Kaspi, & Woods 
2002), further strengthen its classification as an AXP. 

XTE J1810— 197 has a 5.54 s period, and an unsteady 
spindown characterized by a P m 10~ n s s _1 period 
derivative measured in the X-rays (Ibrahim et al. 2004; 
Gotthelf & Halpern 2005). The detection of radio emis- 
sion from the source, first ever for an AXP (Camilo et 
al. 2006), allowed for a more closely-spaced monitoring 
of its period and resulted in the measurement of a larger 
range of spindown rates (Camilo et al. 2007). The range 
of Z?dip corresponding to the observed period derivatives 
are used in Figure 5. 

In an earlier investigation, we reported on the physical 
evolution of this source during its decline from outburst 
(Giiver et al. 2007). Here, we focus on the magnetic field 
strength of XTE J1810-197 using the observation with 
the highest number of counts. 

DISCUSSION 

Figure 5 shows the magnetic field strengths obtained 
for the 4 AXPs by fitting their X-ray spectra with the 
STEMS model against the dipole field strengths inferred 
from the spindown of these sources according to the 
dipole spindown formula. The error bars correspond to 



TABLE 2. Spectral Results of STEMS Model for 4 AXPs 



Source 


Magnetic Field 


Surface Temperature 


T 


P 


*v(d.o.f.) 


4U 0142+61 
IE 1048.1-5937 
IRXS J1708-4009 
XTE J1810-197 


4.60+0.07 
2.26+0.05 
3.96+0.17 
2.68+0.06 


0.31+0.01 
0.37+0.01 
0.35+0.01 
0.31+0.01 


3.54+0.14 
3.91+0.52 
5.26+0.36 
2.36+0.36 


0.43 ±0.01 
0.22 ±0.02 
0.48 ±0.01 
0.25 ±0.02 


0.931 (462) 
0.952 (611) 
1.050(1187) 
1.07 (732) 
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FIGURE 1. STEMS model fit to the X-ray spectrum of 
4U 0142+61. 
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FIGURE 2. STEMS model fit to the X-ray spectrum of 
IE 1048.1-5937. 
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FIGURE 3. STEMS model fit to the X-ray spectrum of 
RXS J1708-4009. 
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FIGURE 4. STEMS model fit to the X-ray spectrum of 
XTE J1810— 197. 



2 — a uncertainty in the values of the spectroscopic mag- 
netic field strength, while error bars on the dipole spin- 
down field reflect the range obtained from the variable 
period derivatives seen in some sources. 

In the cases of IE 1048.1-5937, IRXS J 1708 -4009, 
and XTE J1810— 197, we find a very good agreement, at 
an unexpected level, between the spectroscopically mea- 
sured magnetic field strength and that obtained from their 
spindown. In all three cases, the two values are consistent 
within the formal and expected systematic uncertainties. 
For the case of 4U 0142+61, the spectroscopically mea- 
sured surface magnetic field is a factor of 3 larger than 



the spindown field. This may be due to the simplified 
assumptions in either of the two measurements of the 
magnetic field. Alternatively, this might be an indication 
of multipole magnetic field components on the neutron 
star surface which contribute negligibly to the spindown 
torques. 

The good agreement between the theoretical models 
and the spectral data provide us with a new tool with 
which to understand the physical conditions of mag- 
netar surfaces and magnetospheres. At the same time, 
our spectroscopic measurements of the magnetic field 
strengths offer independent confirmation for the magne- 
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Surface Magnetic Field (x 10 14 G) 

FIGURE 5. The comparison of the spectroscopically mea- 
sured magnetic field strengths of five AXPs to the dipole fields 
inferred from the spindown properties of these sources. The 
error bars in Bjjp represent the range of measured spindown 
rates for each source, while the error bars in the spectroscopic 
magnetic field strength represent 2 — 0" statistical uncertainties. 



tar nature of AXPs. 
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